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Abstract 

Drell-Yan production in the forward region at the Large Hadron Collider is sensitive to multiple 
radiation of QCD partons not collinearly ordered, emitted over large rapidity intervals. We propose 
a method to take account of these radiative contributions via a factorization formula which depends 
on the unintegrated, or transverse momentum dependent, splitting function associated with the 
evolution of the initial-state sea quark distribution. We analyze this formula numerically, and point 
out kinematic effects from the initial-state transverse momentum on the vector boson spectrum. 

> 

in 

r-> 

in 
o 

<N 



1 



I. INTRODUCTION 



Many aspects of the experimental program at the Large Hadron Collider (LHC) de- 
pend on the analysis of processes at large momentum transfers containing multiple hard 
scales. As such, they will be influenced by improved formulations of factorization in QCD 
at unintegrated level pQ, which serve in multi-scale processes to control perturbative large 
logarithms to higher orders of perturbation theory and to describe appropriately nonper- 
turbative physics in the initial and final states of the collision. These formulations involve 
transverse-momentum dependent (TMD), or unintegrated, parton decay and parton density 
functions (pdfs) [2]. 

A broad class of such multiple-scale events is given by small-x processes. These processes 
are one of the main sources of final states in the central region at the LHC [3]. Besides, they 
are responsible for the sizeable rates of forward large-p^ jet production at the LHC, giving 
rise to new phenomenology in the LHC forward region compared to previous collider exper- 
iments, e.g. forward jet physics [HE]. Small- a; TMD factorization and pdfs (see for instance 
recent reviews in [BJ U\) serve for the perturbative resummation of high-energy logarithmic 
corrections and for the development of parton shower algorithms capable of incorporating 
multi-gluon coherence including the effects of large-angle, noncollinear emission. 

Most of TMD computational tools relevant to small-x and forward physics at the LHC 
have so far been developed within a quenched approximation in which only gluon and valence 
quark contributions are taken into account at TMD level [SJEIIS]- While this gives a sensible 
approximation, based on the dominance of spin-1 exchange processes at high energies, to 
the asymptotic behavior of production processes coupled to gluons such as heavy flavor and 
scalar production, it is mandatory to go beyond this approximation to include preasymptotic 
effects, and to treat final states associated with quark-initiated processes such as Drell-Yan 
production. In this work we take first steps to address this issue by including sea quark 
contributions and examine forward Drell-Yan production. 

Drell-Yan processes at the LHC are instrumental in precision electroweak measurements, 
in luminosity monitoring and pdf determinations, and in new physics searches [3] . Studies of 
Drell-Yan at TMD level have recently been performed in the framework of the soft-collinear 
effective theory [9HTT]. Results on high-energy logarithmic corrections to Drell-Yan have 
been obtained in the qg* channel [121 H3] and in the associated production channel Z/W+ 
heavy quarks [13] . On the other hand, early attempts [T5HT8] to treat the unintegrated pdf 
evolution beyond the quenched approximation include quarks via splitting probabilities to 
lowest order of perturbation theory, neglecting any transverse momentum dependence in 
the branching. In [IH] k^-dependent kinematic corrections are included, while the splitting 
kernels are still taken in lowest order. Also, a program to perform shower Monte Carlo 
evolution at unintegrated level has recently been proposed [20] based on the expansion [21] 
in two-particle irreducible (2PI) kernels. This program is formulated at the next-to-leading 
logarithmic order. However, this does not include small- x logarithmic effects which are 
present beyond NLO in flavor singlet distributions. 

The approach of the present paper is based on the high-energy form of the (off-shell) 
TMD quark Green function introduced in [22]. This is obtained by generalizing to finite 
transverse momenta, in the high-energy region, the 2PI expansion [21J. The main point 
is to construct unintegrated sea-quark distributions incorporating the effects of the TMD 
gluon-to-quark splitting kernel |22j, which contains all single-logarithmic small- x corrections 
to sea quark evolution for any order of perturbation theory. In order to relate this parton 
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splitting kernel to forward vector boson production, we analyze the flavor exchange process 
at high energy according to the "reggeized quark" calculus [23] [23] . This extends the effective 
action formalism [23], currently explored at NLO |26j, to amplitudes with quark exchange in 
terms of effective degrees of freedom, the so-called reggeized quarks [23 I2B] • The use of the 
effective vertices [231 [23] ensures gauge invariance of the coefficients relevant to perform the 
high-energy factorization [22} [29] for vector boson production, despite the off-shell parton. 

We then examine kinematic effects related to the difference between the virtuality of the 
exchanged parton and its transverse part which, although formally subdominant both in the 
collinear and in the high-energy expansions, can nevertheless be numerically non-negligible. 
These kinematic effects are similar to those observed in the study of [19]. We perform a 
numerical study of the kinematic terms at the level of the partonic Z cross section, and 
compare this with the qg* — > Zq result [12]. 

The formulation proposed in this work can be implemented in a parton shower Monte 
Carlo event generator. Details of the formulation and shower implementation will be de- 
scribed elsewhere [30] . The outline of the present paper is the following. In Sec. 2 we apply 
the reggeized quark calculus to examine Z boson production and reobtain the TMD quark 
splitting function. In Sec. 3 we investigate kinematic effects in the factorization formula and 
study this numerically. We give conclusions in Sec. 4. 



II. REGGEIZED QUARK AND UNINTEGRATED QUARK DENSITY 



In this section we relate the Z boson production cross section to the unintegrated quark 
density defined at high energy via the transverse momentum dependent kernel [22]- For 
definiteness we consider factorization of the g*q — » Zq matrix element into off-shell point- 
like qq* Z coefficient and off-shell gluon-to-quark splitting function, which is of direct 
phenomenological relevance for forward Drell-Yan production. 

Within the reggeized quark formalism of [23] [23] this process can be described in the 
high energy limit by the single effective diagram in fig. [I] The double lines (with fermionic 
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FIG. 1. The g*q —> Zq process within the reggeized quark formalism. Double lines with arrow 
indicate the effective reggeized quark exchange in the t-channel. 



arrows) in the diagram represent the exchange of off-shell (reggeized) quarks, while the black 
dots represent their effective couplings to s-channel quarks and gauge bosons [23]. We use 
the following Sudakov decomposition for the initial and final momenta, 

,2' 



V = X2P2 



k = X1P1 + k 



q = zxxV\ + [x 2 + 
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where pi and p 2 are light-like four-momenta respectively in the plus and minus lightcone 
directions, 2pi ■ p% — s, and C 2 = t. The reggeized quark propagators are purely transverse 
fermion propagators, supplemented with a projector in longitudinal momentum space, 



T i 



2pi ■ P2 ~£ 2 + it 



2pi-P2 -I 2 



IE 



(2) 



The couplings of the off-shell quarks to Z bosons and usual quarks are 



p q,i* 



cos 9 W sin 9w 



r»(j>,£,q) (V f -A fl5 ), 



(3) 



with the "plus" off-shell reggeized vertex given by 



Pi 

Pi ■ q 



e- Pl =o 



(4) 



where the light-cone pi component of the t-channel momentum i is set to zero. Similarly, 
the lower vertex is given by 



tgt a T»_(kJ,p') 



(5) 



with 



P2 

p 2 ■ k 



(6) 



where the light-cone p\ component of the t-channel momentum is set to zero. 

The kinematic conditions I ■ p\ = and I ■ p 2 = in eqs. Q and ^ reflect the strong 
ordering in the minus and plus lightcone momenta j2H |25] and are natural from the point 
of view of high-energy factorization. Using these conditions, the expression which results 
from fig. [T] agrees precisely with the high-energy expansion of the g*q — > Zq matrix element, 
as carried out in [12] along the lines of the high-energy resummation [22J. However, the 
expression resulting from the strong ordering kinematics in eqs. Q and only yields a 
good approximation in the limit of asymptotically large partonic center-of-mass energy. In 
the present analysis, we concentrate on a different aspect. On one hand, we are interested in 
improving the description for center-of-mass energies which are not asymptotically large. On 
the other hand, we are interested in obtaining an off-shell (i.e., /cy-dependent) factorization 
in terms of the transverse-momentum dependent sea-quark distribution defined in [22] from 
the two-particle irreducible expansion. 

To this end, we will relax some of these kinematic approximations. In the rest of this 
section we match the high-energy factorized expression to collinear factorization by relaxing 
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the condition £ ■ p 2 = 0. In Sec. 3 we discuss the effect of relaxing the condition £ ■ pi = 0. 
Instead of eq. ^ we consider the effective vertex with exact kinematics 



rt(k,l,p')=igt a [Y + -™-f 



l-l 
P'2 



with 



-ZX\P\ 



t + q 2 



p 2 • k J ZX\S 

Current conservation is satisfied j2U [25] also in this more general case 

fc M ft(M,p>(z/) = igt a ^u(p) = 0. 



-P2 - q 



(7) 



(8) 



In order to combine the above vertices with gluonic kr factorization [29], we contract the 
above vertex with the transverse gluon momentum, 



rt(k,£, P ') 



fe, 



'k 2 



(9) 



This yields the q*qg*-vertex which will be used in the following. Owing to eq. (|8j), this is 
merely a rewriting from the longitudinal polarization naturally associated with reggeized 
gluons to the transverse polarization, which allows one to perform the matching with the 
collinear factorization [211 E2] • Removing the condition £ ■ p 2 = for the effective vertex, we 
also promote the propagators in eq. ^ to full four-dimensional propagators as follows, 



£ 2 + ie 2pi ■ p 2 



i-f 



2p 2 ■ pi £ 2 + it 



(10) 



As in [2H 122] , to perform the matching in the collinear region we include an upper bound 
on four- momentum exchanged in the i-channel, £ 2 < fi 2 F , with /xp the factorization scale. 

Using the setup described above, we find that the kx factorized qg* — > Zq cross-section 
can be written as 



kT— fact. 



dz J -^®{V 2 F-? 2 ) H^2S,M 2 z ,z,q 2 ) R g * g *(z,q,k), (11) 



where a is the following off-shell continuation of the pointlike qq — > Z matrix element, 

a(x lX2 s, Mf, z, q 2 ) = V2G F M 2 Z (V 2 + A 2 q )^6(z Xl x 2 s - q 2 - M|), (12) 

with Mz the mass of the Z-boson, Gf the Fermi coupling, V q and A q the coefficients of the 
vector and axial couplings of the Z to a quark, and R g * q * (z, q, k) is given by 



R g * q *(z,q,k) 



T R a s 



z(l-z) 



2tt (A 2 + z(l - z)k 2 ) 



A 2 (A ■ k) 2 

+ 4(l-2^)A-fc-4 l ] 



z 1 - z 



k 2 



+ 4z(l - z)k 2 



(13) 



where Tr = 1/2, A = q — zk. Upon azimuthal average, eq. (13) returns the TMD gluon-to- 



quark splitting kernel [22] • We examine azimuthal effects numerically in the next section. By 
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replacing R g * q *(z, q, k) with its azimuthal average and neglecting the order-zk contribution 
to the transverse momentum recoil in eq. (12), we get 



fcr-fact. ^ 
°qg*^Zq - 



where 



dz 



d 2 A 
VA 2 



a(xix 2 s, M§, z, A 2 ) 6 ( - zk 2 - 
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P qg {z,k 2 ,A 2 ) 



P qg (z,k 2 ,A 2 )=T R 



A 2 + z(l - z)te 



2 r 



z) 2 + z 2 



Az 2 (l 



(14) 



(15) 



It was pointed out in [22J that eq. (15) allows one to relate the sea quark distribution, 



with next-to-leading logarithmic accuracy at small x, to the small- x gluon Green's function 
Q (x, k 2 , ii 2 ) obeying the BFKL equation. Although evaluated off-shell, the splitting function 



in eq. (15) is universal [221 IS]- It reduces to the collinear splitting function at lowest order 
for k = 0, but in the high energy limit it factorizes correctly the finite k dependence to all 
orders. By using this result, the unintegrated sea quark distribution is given by 



X. 



\i 2 ' /i 2 



dz 

z 



dk 2 /4- 



A 2 + z(l - zW 



1 a s 
A 2 2n 



-P qg (z,k 2 ,A 2 )g(^k 2 ^ 2 ). (16) 
We thus write the forward Z-boson hadronic cross section in the high energy factorized form 



<j, 



pp- 



>z(s 



Mf) 



J dxi J ~~ J dx 2 cr Q s 



n 2 



A 2 ii\ 



Q 



(?) 



fl 2 



(17) 



with Q sea given in eq. (16) and the (integrated) parton distribution as provided by 



collinear factorization. The scale \i in eqs. (16) and (17) is an initial infra-red scale at 



which the parton distribution functions are defined. The dependence on this scale occurs 
for the TMD distributions as a direct consequence of the matching procedure between high 
energy and collinear factorization. The presence of collinear singularities in the integrations 
over momentum lines connecting the 2PI amplitudes requires the use of a regulator for the 
BFKL gluon Green's function which breaks the scale invariance of the LO BFKL equation. In 



addition to the BFKL gluon density, also the integral over A 2 in eq. (17) is collinear divergent 
in the limit k 2 — > and requires regularization, leading to an additional dependence on 
the scale fi 2 . While in [22] this is achieved through dimensional regularization in d = 



4 + 2e dimensions, a cut-off regularization will be used in future Monte-Carlo realizations 
of eq. (16). The dependence on the scale /x cancels up to second order corrections, i.e. 
order a s (a s In l/x) n which is the accuracy at which eqs. (16) and (17) are valid. 



to 



We will be interested in applications of the sea quark distribution in eq. ( 16 ) to shower 



Monte Carlo generators including the transverse momentum dependence of the parton 
branching, such as [32]. Compared to previous TMD approaches [151 EH I2H] , this distribution 
includes the finite-fc dependence of the kernel (15), responsible for all-order logarithmically 
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enhanced corrections to quark evolution at small x. For these applications it is relevant to 
investigate the size of corrections associated with the high-energy kinematic approximations 
discussed earlier. We turn to this in the next section, and study the numerical effect of the 
kinematic contributions on the off-shell Z-boson production cross section. 



III. KINEMATIC EFFECTS IN THE PARTONIC Z CROSS SECTION 

In this section we discuss that the accuracy of the factorized expression obtained above 
can be increased by going beyond the strong ordering approximation in eqs. ([3]),(|4]). We will 
write the convolution integral in terms of the virtuality of the four-momentum exchanged in 
the t-channel, \t\ = —£ 2 , rather than the transverse momentum. While the two formulations 
coincide in the asymptotic high-energy limit z — > 0, they differ by inclusion of finite- z terms. 
We examine numerically the role of these terms. 

The formulation in\t\ arises naturally if we remove the constraint £-p% = in the effective 
vertex of eq. Q. The kinematically improved version of this vertex is 

tW,q)= (V+/-— V (18) 

V Pi-qJ 



Gauge invariance can be verified similarly to eqs. ([8]), ([9]). One then recovers eq. (11) with a 
replaced by 



& = V2G F M 2 z (V q 2 + A*)^-6(z Xl x 2 s + t - Mf), (19) 



where 



-t = -t = + zk 2 . (20) 

1 — z 

The angular average can now be performed exactly, and yields 



1 oo 



, r f dA 2 tn { A 2 + z(l-z)k 2 \ a s . ,„ a9N 
< 9 ^z q = J dz J ° (tf? 1 _ z ) ^ P « 9 {*> k ^ A2 ) • (21) 



Unlike eq. (14) the above expression no longer decouples in transverse momentum space. 



However, by changing variable to \t\ in the convolution integral, we get 

i 

vl g *^z q ( x i x 2S,k 2 ,Mz) = J dz J d ^\ k 2 6 (/4 - N) cr{xiX 2 s,Ml,z,\t\) 



x Q(\t\-zk 2 ) ^- Pqg ( z ,k 2 ,(l-z)(\t\-zk 2 )). (22) 

This result shares some features with the results in [17] and (THj. In particular, as in [19] 
it takes into account subleading corrections to strong ordering by keeping track of the ex- 
act kinematic relation between virtuality and transverse momentum. It was shown in the 
numerical NLO-DGLAP analysis of [19] that these kinematic contributions provide a large 



fraction of the full corrections at NLO. On the other hand, unlike jTTJ [19] eq. (22) contains 
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the dynamical effects of transverse momentum dependent terms in the gluon-to-quark split- 
ting kernel. These become potentially important in the forward production kinematics, as 
they give rise to logarithmically enhanced terms for small x at higher orders in a s . 



While the general structure of the unintegrated sea quark density eq. (16) remains un 



changed, the off-shellness is now expressed in terms of the absolute value of the four mo- 
mentum square exchanged in the t-channel, 



Q? a I x 



' 9 ' 9 



x 



In [22l|3l] the low-fe 2 behavior of eqs. (|16|) , (|23|) is analyzed using dimensional regularization. 
For the parton-shower applications [SU] an infrared cut-off is applied. In the numerical 
study that follows, we examine the effect of the TMD sea quark distribution on the Z-boson 
partonic cross section. For this purpose we will work at fixed k 2 . 



We perform a numerical comparison of eq. (14) and eq. (22), corresponding respectively 



to the kinematic effects of going beyond the plus-momentum strong ordering and the minus- 
momentum strong ordering in the high-energy factorized formula. We compare this with the 
full qg* — > qZ matrix element result. We consider also the result of the collinear approxima- 
tion obtained by setting k 2 
momentum variable 



in eq. (14). We examine the dependence on the transverse 
zk. 



A = q — zk. Here q is the vector boson transverse momentum; z is 
the plus-momentum fraction transferred to the vector boson from the incoming gluon that 
results from the small-x initial-state shower; k is the transverse momentum carried by the 
partons emitted in this shower (i.e., partons radiated in addition to the leading quark against 
which the vector boson recoils). In fig. [2] and fig. [3] we examine respectively the region where 
|A| is small compared to the Z boson mass Mz and the region where |A| is on the order of 
the Z boson mass or larger. 

For small |A| we find that the differences between eq. (14) and eq. (22) are numerically 



small, and that both expressions are close to the full result; as |A| increases, we find that the 



deviations due to the kinematic contributions by which eq. ( 14 ) and eq. ( 22 ) differ become 
non-negligible, and that eq. (22) gives a better approximation to the full result. Fig. [2] 
illustrates the small |A| region. In fig. |2](a) we plot the differential cross section da/dA. 2 for 
fixed values of k 2 and X1X2S. The solid, dashed and dot-dashed curves overlap on the scale 
of this plot, while the dotted curve, corresponding to the collinear approximation, deviates 
from them when |A| < |fc|. In fig. [2^b) we zoom in on the small relative deviations of 
eq. (14) and eq. (22) from the full result for the differential cross section da/dA. 2 . We see 
that throughout the range shown in fig. |27b) the deviations are at most of the order of few 



percent. In this region the factorized expressions based on reggeized quark graphs are close 
to the full result, independently of the different kinematic approximations on the longitudinal 
momentum orderings. What dominates this region are the transverse momentum dependent 



terms, which are kept correctly by both eq. (14) and eq. (22). On the other hand, we see in 
fig. |2|(a ) that in the small |A| region the collinear expression obtained by setting k 2 = in 



eq. (14) is far from the full result, due to missing TMD corrections to the kernel (15). 

As |A| increases the kinematic effects from the plus-momentum ordering and minus- 
momentum ordering become significant. Fig. [3] illustrates these effects for |A| on the order 
of the Z boson mass Mz, or larger. Note that in this region the collinear k 2 = expression 
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FIG. 2. (a): A 2 dependence of the differential cross section da/dA 2 for small |A|: (solid) full; 
(dashed) no plus-momentum ordering; (dot-dashed) no plus-momentum and minus-momentum 
ordering; (dotted) collinear approximation. All but the last curve overlap in this region. We set 
x±X2S = 2.5M|, k 2 = 2 GeV 2 . (b): Relative deviations in the differential cross section da/dA 2 : 
(dashed) no plus-momentum ordering; (dot-dashed) no plus-momentum and minus-momentum 
ordering. 



gives a good approximation to eq. (14). The dashed and dotted curves nearly overlap on 
the scale of the plots in fig. [3j However, the full result contains significant corrections with 
respect to these curves, corresponding to terms that are subdominant both in the collinear 
and in the high-energy expansions. We see that a non-negligible part of these corrections 
can be taken into account by including kinematic contributions from terms beyond the 
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do- 1 

l(T 12 x [ ] 

dA 2 GeV 4 
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(b) 

FIG. 3. A 2 dependence of the differential cross section da/dA 2 in the large |A| region: (solid) 
full; (dashed) no plus-momentum ordering; (dot-dashed) no plus-momentum and minus-momentum 
ordering; (dotted) collinear approximation. We set k 2 = 2 GeV 2 ; (a) X1X2S = 2.5M|, (b) X1X2S = 
12.5M 2 . 



strong ordering in the lightcone minus momenta, as eq. (22) does 
effects observed in [15] . As a result eq 
matrix element 



This is analogous to 
(22) provides a closer approximation to the full 



It can be verified numerically that the corrections die out in the limit of 
asymptotically large s. This analysis is however significant for applications to finite collider 
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energies. 

We thus observe that taking into account TMD terms in the parton splitting and the 
reggeized quark coefficient enables one to extend the description of vector boson production 
in terms of factorized quark distributions from the collinear region into the region where | A| 
is of order \k\ and transverse momenta are no longer strongly ordered. This is a region where 
collinear approximations are seen not to be sufficient, and it is relevant to forward Drell-Yan. 
We will use the unintegrated quark distribution thus defined in forthcoming studies [30J . As 
a result of TMD terms, on the other hand, the initial state kinematics becomes rather more 
complex. Corrections to longitudinal momentum ordering in the effective vertices are seen 
to be numerically significant when |A| is of order the vector boson mass. To this end it 
is useful to employ the formulation given in this section which includes the exact one-loop 
kinematics. 

IV. CONCLUSION 

Drell-Yan production processes are instrumental in hadron collider experiments both for 
hadronic physics studies and for new physics searches. At the LHC a new kinematic region 
opens up for forward Drell-Yan, in which the QCD treatment of the production process and 
associated final states acquires new features due to multi-parton radiation over long rapidity 
intervals. The emission of soft gluons that are not collinearly ordered becomes relevant, and 
effects of color coherence set in associated with the region of small longitudinal momentum 
fractions x. Taking this into account calls for improved formulations of factorization in QCD 
at unintegrated, or transverse momentum dependent, level. 

Fully general TMD factorization formulas are still lacking [U |2] , due to the difficulty in 
disentangling systematically soft and collinear gluon correlations between initial and final 
states. This is illustrated, for instance, by back-to-back di-hadron and di-jet hadroproduc- 
tion [331 El]- In the case of small x, however, TMD factorization results exist [23122] based 
on dominance of single gluon helicity at high energy. These results can be used to define 
gauge- invariantly the TMD gluon distribution. See recent works in [E1E51EE]- In this paper 
we have used the method [22] to extend the above treatment to the TMD sea-quark distri- 
bution. This is the dominant channel coupling to forward Drell-Yan production. We have 
related the TMD splitting kernel [22] to the forward Z-boson cross section. The reggeized 
quark formalism [231 121] is used to treat the off-shell quark coefficient and ensure gauge in- 
variance. We have analyzed numerically effects due to the transverse momentum kinematics 
in the relationship between the Z boson spectrum and the TMD sea-quark distribution. 

Phenomenological applications of the reggeized quark formalism were investigated in [2H] 
using the approach [161 02] to TMD pdfs. Compared to these studies, the main feature of our 
approach is that it includes, besides the lowest-order quark splitting function, the full series 
of finite-fc_i_ terms in the TMD kernel. These terms control the perturbative summation of 
small-x logarithmic corrections to flavor-singlet observables to all orders in a s . 

Future extensions of the results above involve several directions. One concerns large- 
x corrections, likely to be important for Drell-Yan phenomenology. See [HI [3TH3U] for 
discussion of x — > 1 issues in TMD quark distributions. Another direction concerns nonlinear 
effects from high parton densities in the small x region. Recent works [iOll^T] study multiple 
scattering contributions in dense targets and nuclei. Techniques such as [12] have been 
proposed to incorporate the treatment of multiple-gluon rescattering graphs at small x 
starting from the operator matrix elements [H |2] for parton distributions. They may be 
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helpful for extensions to the high density region that are aimed at retaining accuracy also 
in the treatment of contributions from high px processes, e.g. Drell-Yan on nuclei [43j. 

The results in this paper can be implemented in a parton shower Monte Carlo generator 
including transverse-momentum dependent branching, such as [32] . Work along these lines is 
in progress [30] . This will allow one to address predictions for exclusive final-state observables 
associated with forward Drell-Yan production, for which intense experimental activity is 
forthcoming [44J. 
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